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Abstract 

The critical behavior of frustrated Josephson-junction arrays at / — 1/2 flux quantum per pla- 
quette is considered. Results from Monte Carlo simulations and transfer matrix computations support 
the identification of the critical behavior of the square and triangular classical arrays and the one- 
dimensional quantum ladder with the universality class of the XY-Ising model. In the quantum 
ladder, the transition can happen either as a simultaneous ordering of the Z2 and (7(1) order param- 
eters or in two separate stages, depending on the ratio between interchain and intrachain Josephson 
couplings. For the classical arrays, weak random plaquette disorder acts like a random field and 
positional disorder as random bonds on the Z2 variables. Increasing positional disorder decouples the 
Z2 and (7(1) variables leading to the same critical behavior as for integer /. 
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1 Introduction 



For an ordered two-dimensional array of Josephson junctions at rational values of the flux quantum per 
plaquette / — ^/^o, where $ is the magnetic flux through a plaquette and $o the flux quantum, the 
ground state is a pinned vortex lattice commensurate with the underlying periodic pinning potential 
leading to discrete symmetries in addition to the U{1) symmetry of the superconducting order parameter 
[|]. When quantum fluctuations due to charging effects can be ignored, yielding the so-called classical 
array, the critical properties can be described by a frustrated XY spin model in two dimensions with 
frustration parameter /. Inclusion of charging effects adds an extra dimension, the time direction, but one- 
dimensional arrays in the form of ladders can still be described by two-dimensional classical spin models. 
Of particular interest is the case when / = 1/2 which has been intensively studied both experimentally 
and theoretically because of the interplay between continuous U{1) and discrete Z2 order parameters. 
The Z2 Ising symmetry can be associated with an antiferromagnetic arrangement of plaquette chiralities, 
X = ±1, measuring the directions of the circulating currents in each plaquette ^ |j. 

It is quite natural to expect for a Josephson-j unction array ai f — 1/2, that critical behavior resulting 
from the interplay between vortices and Ising degrees of freedom, can be described by some form of coupled 
XY and Ising models. Such models, containing both continuous U{1) and discrete Z2 symmetries, can 
give rise to interesting and sometimes unusual critical behavior Q]. Recently, there has been a lot 
of interest in these systems, not only in relation to Josephson-junction arrays but also in surface phase 
transitions where the interplay between roughening and reconstruction of the surface is important ^ ^ . 
A minimal model which is general enough and consistent with the basic symmetries can be defined by 
the XY-lsing model [|, |, § 

^ = - II [(^ + Ba,a,) cos{e, -6,)+ Ca,a,\ (1) 
<ij> 

where a = ±1 is an Ising spin and 9 = [0, 27r] is the phase of a two-component unit vector [XY spin). The 
Ising and XY variables are coupled by their energy densities in the same way as two independent systems 
of Ising spins are coupled in the Ashkin- Teller model. The model can also be regarded as the infinite 
coupling limit, /i — > 00, of two XY models |T^, ^ |l2| coupled by a term of the form —ft- cos 2(6*1 ^6*2)- 
In the j4i?-plane, the model defined by Eq.^) has a rich phase diagram that depends strongly on the 
value of C. The model with B is relevant for an anisotropic frustrated XY model 12. 1^ which 



can, in principle, be physically realized as a modulated Josephson-junction array at / = 1/2, where the 
Josephson coupling in every other column of junctions is different by a constant factor from the others. 

For the isotropic square and triangular classical arrays and for the ladder quantum array [l5[| the 
relevant subspace of model (|l|) is the symmetric one, A ^ B, 

<ij> 

This case has a special symmetry since XY variables are not coupled across an Ising domain wall where 
CifTj -f 1 0. In the original frustrated XY model, the phases across the chiral domain walls are actually 
coupled in the ground state This coupling and additional terms allowed by symmetry have always 
been assumed to be irrelevant at criticality in the XY-ls\ng model description. Recently, Knops et al [|6j 
have addressed this question using transfer matrix calculations and provided further justification for the 
relation between the XY-lsvcig model and the frustrated XY model on a square lattice by showing that 
at the transition the phase coupling across domain walls vanishes for increasing system sizes. 

The phase diagram of the XF- Ising model given by Eq. (||) as inferred by Migdal-Kadanoff renormal- 
ization and Monte Carlo simulations ^, || consists of three branches which meet at multi-critical point 
P, as indicated in Fig. ^ One of the branches, PT, corresponds to single transitions with simultaneous 
loss of XY and Ising order, and the other two to separate Kosterlitz-Thouless and Ising transitions. 
The line of single transitions becomes first order at the tricritical point T. The critical behavior of a 
Josephson-junction array at / = 1/2 as a function of temperature corresponds to a particular cut through 
this phase diagram and the single or double character of the transition depends on the relative location of 
this path to the branch point P. In Monte Carlo simulations the critical line PT in the phase diagram 
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appeared to be non- universal as the critical exponents associated with the Z2 order parameter were found 
to vary systematically along this line. More recent work based on transfer matrix calculations |]l7| found 
no clear evidence for variation of these exponents and suggested that the apparent non-universality may 
be due to slow convergence and strong corrections to scaling. However, it is clear that the values of the 
critical exponents estimated by finite-size scaling of large systems are inconsistent with pure Ising critical 
behavior, = 1, 77 = 1/4, as can be seen from Table 0. 

Table 1: Critical exponents v and r\ for the Z2 order parameter and central charge c of the Xy-Ising 
model on the critical line PT (Fig. |l|) at A = 1, C = —0.2885. Results are from Monte Carlo simulations 
(MC) and Monte Carlo transfer matrix calculations (MCTM). The value of c is the estimate for the 
largest system without extrapolation 







V 




c 


MC 


Ref. PI 


0.84 


0.31 




MCTM 


Rcf. H 


0.79 


0.40 


- 1.60 



The central charge c, which provides additional information on the nature of this critical line of 
simultaneous Ising and XY transitions was also evaluated by transfer matrix calculations on infinite 
strips of widths as large as 30 lattice spacings |l^. A value, c = 1 -I- 1/2, would be expected if the critical 
behavior is the result of a superposition of critical Ising [c— 1/2) and Gaussian ic—V) models |l^. The 
numerical estimate for c however is usually higher than this value, as indicated in Table ^ but effective 
values of c obtained for increasing system sizes decrease significantly. Extrapolation assuming power-law 
corrections |17 yield results not inconsistent with c = 3/2. These results suggest that the critical line PT 
is controlled by a non-trivial fixed point with new critical exponents. In this case the phase transition 
in the Josephson-junction array system at / = 1/2 would be in a new universality class, if this system 
indeed corresponds to a cut through the XF-Ising model that intersects the critical line of simultaneous 
XY and Ising transitions. 



2 Classical arrays: the frustrated XY model 

There have been several numerical studies of the frustrated XY model on triangular and square lattices, 
but so far no definitive conclusion has been reached as to the nature of the transition. These are the 
simplest models for a two-dimensional array at / = 1/2 when capacitive effects can be neglected [0 
and the critical behavior describes the superconducting to normal transition in these systems. Different 
numerical techniques, including finite-size analysis, seem to lead to conflicting results which suggest 



either a single |19|, |2C|, ^ ^ |16| or double transitions ||23|, |2J, |25[ ^ . It is clear however, that if there 
is no single transition then there are two transitions which occur at temperatures very close to each 
other and in a generalized version of the model they could join depending on the effective parameter 
controlling the XY and Ising excitations. In fact, a Coulomb-gas representation of the frustrated XY 
model with an additional coupling between nearest-neighbor vortices suggests a mechanism that reconciles 
the apparently contradictory results: this model has a phase diagram |^ with similar structure as in Fig. 
for the XF- Ising model with the occurrence of double or single transitions depending on this additional 
coupling. Possibly a Josephson-junction array at / = 1/2 corresponds to a cut through the phase diagram 
of the XF- Ising model close to the branch point P . If the critical behavior associated with the C/(l) and 
Z2 symmetries are determined separately, the single or double nature of the transition is likely to remain 
unresolved on purely numerical grounds. Computing the critical point of a Kosterlitz-Thouless transition 
even in the simplest case is problematic and requires scaling forms that lead to highly unstable fits. As 
a consequence, the error bars of the critical point make resolving two critical points impossible in many 
cases. 

If the frustrated XY models and the coupled XYAsmg model are in fact in the same universality 
class then in order to verify the single nature of the transition it is sufficient to study the Z2 degrees of 
freedom M, EOl. If the critical exponents are inconsistent with pure Ising values, the transition cannot 
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correspond to the Ising branch of a double transition as in Fig. |1[ This point of view has the advantage 
of not requiring a precise study of the behavior of the XY variables which for a model of this nature 
cannot be interpreted using what is known from the standard XY model. The Z2 critical exponents of 
the XY-ls'mg model on the critical line have been compared with those of the frustrated XY model on 
square and triangular lattices using the same numerical methods po|] . In particular, the estimate of the 
thermal exponent, ly, was obtained by finite-size scaling of Monte Carlo data ||20|| which is insensitive 
to the estimate of the critical temperature, T^. In most of the other numerical work, the exponent v 
is determined in an indirect way with several fitting parameters while this evaluation involved only a 
one-parameter fit to the data. The results, f = 0.83, 77 — 0.28 and = 0.85, rj — 0.31 for the triangular 
and square lattice, respectively, are significantly different from the pure Ising values and thus favor a 
single-transition scenario. They are also in fair agreement with the critical exponents in Table ^ for 
the XY-lsing model on the critical line. These were the first results which clearly indicated that the 
chiral exponents deviate from pure Ising critical behavior. Independent estimates using finite-size scaling 
of correlation functions ||2^, Monte Carlo transfer matrix calculations ||2l[, simulations of the lattice 
coulomb gas ||2^ and exact transfer matrix calculations seems to agree on the estimate of v and the 
deviation of the chiral exponents from the pure Ising behavior. 

The deviation from the pure Ising exponents for the square lattice has recently been questioned p6| on 
the basis of a possible failure of the finite-size scaling used in previous estimates. Under the assumption 
that the transition takes place in two stages very close in temperature, another divergent length scale 
should be included in the scaling analysis which could be the reason for the non-Ising exponents reported 
so far. The result for the temperature dependence of the Z2 correlation length appear to be consistent 
with the pure Ising exponent v = I but we note that this calculation was performed in a temperature 
range near Tc which is significantly greater than the one where deviations from the Ising exponent have 
been found |2^, |l^. Due to the nearby Ising and XY transitions as in Fig. ^ and crossover effects, a 
possible non-trivial fixed point on the single line can only be reached if calculations are carried out close 
enough to the critical line which requires both large systems and temperatures very close to Tc. Other 
calculations of the Z2 correlation length but closer to Tc are consistent with deviations from the pure 
Ising behavior . 

Table 1^ summarizes recent results obtained for the chiral order parameter of the frustrated XY model. 

Table 2: Recent estimates for the chiral critical exponents v, rj and central charge c of the frustrated XY 
model on a square and triangular(T) lattice. Results considered consistent with v = 1 arc indicated by 
~ 1 



Rcf. jl9| 

Rcf. 1201 

Ref. 1201 

Ref. 1281 

Ref. 1221 

Ref. 1211 

Ref. 1251 

Ref. 116 1 

Ref. 1261 



- 1 
0.82 
0.83 

- 1 
0.86 
0.80 
0.84 
0.77 

- 1 



0.40 
0.31 
0.28 (T) 
0.26 
0.22 
0.38 
0.26 
0.28 



1.66 



1.61 



1.55 



3 Effects of disorder 

In real arrays, disorder is always present due to fabrication processes but it can also be deliberately 
introduced in order to study its effects [ p9| . Two kinds of disorder may be present: random plaquette 
areas and positional disorder |3l] . These types of disorder are irrelevant in zero applied magnetic 
field but become important and have interesting effects at higher fields. Randomness of plaquette areas 
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is realized when the superconducting grains have random sizes but the magnetic field is sufficiently low 
to permit a partial Meissner effect in individuals grain so that the fluxes can be regarded as independent 
random variables. This random-flux type of disorder has been argued to lead to the destruction of 
superconductivity and to a finite correlation length [Eol, ^ ~ which depends on the field H, but 



vortex pinning [S^, 33 by disorder is likely to dominate the superconducting behavior. If only weak 
positional disorder is present, phase coherence is possible within a range of temperature < T < T"*" 
for sufficiently low fields or disorder |30[ For an average value of flux quanta per plaquette, /o, 

the maximum field or disorder for which quasi-long-range phase coherence is possible is bounded by 
/oA < 1/V327r where the probability distribution of the grain positions is P{u) oc exp(— m^/2A^) and u 
is the displacement of a grain from its average position. The analysis also shows that the behavior of the 
superconducting transition at has the same features as the Kosterlitz-Thouless transition but with 
a non-universal jump of the helicity modulus. The disorder-induced reentrant phase transition which is 
predicted to occur at the lower temperature, T^, does not take into account vortex pinning by disorder. 
Since the dynamics is expected to be very slow, experimental or numerical study of the this transition is 
difficult. 

Experiments on proximity-coupled arrays with deliberate positional disorder and numerical sim- 
ulations [^5[ |3^ for integer fo appear to be consistent with the existence of a critical field but find no 
evidence for reentrance. A recent renormalization group treatment | |3^ reaches the conclusion that the 
transition to the normal phase is not reentrant but vanishes at the same critical disorder. 

For fo — n- + 1/2, the ordered array has a Z2 symmetry in addition to the U{1) symmetry. With 
a Coulomb-gas representation it has been shown that random plaquette areas act like random fields on 
the Z2 order parameter and positional disorder acts like random bonds |30[ |. Consequently, random 
plaquette areas should destroy chiral order and phase coherence for any amount of disorder. This result 
seems to agree with energy balance arguments for the chiral ground state |38j. Positional disorder, 
however, is expected to affect the chiral order less dramatically and have different effects on the Z2 
and U{1) excitations. Increasing the amount of disorder should reduce the XY transition temperature 
to below that of the Ising if the Z2 and U{1) degrees of freedom are treated separately. Eventually, 
increasing disorder should lead to a double transition scenario where phase coherence is lost within the 
chiral ordered phase where domain walls and corner charges should have no effect. This leads to the same 
critical field as for integer fo. In the experiments ||29t] on proximity-coupled arrays, the critical field fc{q) 
appears to be independent of q where fo = p/q which is consistent with the proposed scenario |^0|| if it 
is realized for all rational values of fo- Monte Carlo simulations of the frustrated XY model on a square 
lattice finds evidence for the splitting into separated XY and Ising transitions |^ as predicted but no 
evidence for re-entrance [B9l. 



4 Quantum ladder 

A one-dimensional ladder of Josephson junctions, at zero temperature, undergoes a superconductor to 
insulator transition as a function of charging energy Ec due to capacitive effects. The transition is 
strongly affected by the magnetic field |4^, ^ through the value of /. As for two-dimensional arrays, 
the universality class of this transition is a problem of great interest specially in relation to experiments 
[pl 



[15 



^ and theoretical predictions of universal properties 45 . At / = 1/2, the effective Hamiltonian 
l[| describing ffuctuations from a commensurate phase of pinned vortices is expected to be in the 
same universality class of the XY-lsing model in the parameter space A = B. The location of the 
cut through the phase diagram in Fig. ^ depends on the ratio E^/Ey between the interchain E^ and 
intrachain Ey Josephson couplings of the ladder. Since, the superconducting-to-insulator transition at 
T = is to be identified with the loss of phase coherence, this transition in the decoupled region is 
in the Kosterlitz-Thouless universality class, but in the single transition region it is possibly in a new 
universality class, the same as for a two-dimensional classical array at finite temperature. 

The critical behavior of the quantum ladder at / = 1/2 has been studied using a Monte Carlo transfer 
matrix applied to the path-integral representation of the model and a finite-size scaling analysis ||l5|] . In 
this formulation, the one-dimensional quantum model is mapped into a two-dimensional classical model 
with an extra dimension, the time direction. The parameter a = y^Ey/Ec plays the role of an inverse 
temperature in the classical model, where E^ is charging energy. Near the critical point ac, chiral order 
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is destroyed by quantum fluctuations with an energy gap vanishing as jo; — ad''. At the critical point, the 
correlation function decay as a power with an exponent rj. The results for the critical coupling etc and 
critical exponents i' and 77 for two different values of the ratio E^/Ey are indicated in Table ^ For equal 
interchain and intrachain couplings, E^ — Ey, the results for the critical exponents differ significantly 
from the pure two-dimensional Ising model and are in fair agreement with the corresponding values at 
the single line of the XY-lsing model in Table ^ For E^/Ey — 3, the results are in reasonable agreement 
with pure Ising values. In the phase diagram of the XK-Ising model, this corresponds to a path in the 
decoupled region where XY and Ising transitions can take place at different points. In fact, for this 
case, the values obtained for the critical coupling associated with chiral order ai = 1.16, and the one 
associated with the universal jump in the helicity modulus, axY — 1.29, clearly indicate that there are 
two transitions. Somewhere in between E^/Ey — 1 and 3 there should be a bifurcation point where 
a single transition with simultaneous loss of phase coherence and chiral order decouples into separate 
transitions, the chiral transition being located in the insulating phase. 

Table 3: Chiral critical exponents (i^, 77), for the quantum ladder 



Ex/ Ey 


h' 


V 


1 


0.81(4) 


0.47(4) 


3 


1.05(6) 


0.27(3) 



5 Conclusions 

Numerical results from Monte Carlo simulations and transfer matrix computations support the identifi- 
cation of the critical behavior of Josephson-junctions arrays at / = 1/2 flux quantum per plaquette with 
the universality class of the XY-Ising model. For classical arrays without charging effects, in the form of 
square and triangular lattices, the results are consistent with a transition corresponding to a cut through 
the XY-Ising phase diagram intersecting the critical line of simultaneous XY and Ising transitions. Some 
conflicting results regarding the single or double nature of the transition may be explained as resulting 
from the proximity of the cut to the branch point in the phase diagram. For the quantum ladder, the 
ratio between interchain and intrachain Josephson couplings can be used to tune the transition near the 
branch point resulting in a simultaneous or double transition depending on this ratio. For classical arrays, 
weak random plaquette disorder acts like a random field and positional disorder as random bonds on the 
Z2 order parameter. Positional disorder tend to decouple the XY and Ising variables leading to the same 
critical disorder as for integer /. 
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Figure 1: Phase diagram of the XY-lshig model obtained by Monte Carlo simulations jo). Solid and 
dotted lines indicate continuous and first-order transitions, respectively. The precise locations of P and 
T are uncertain. 
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